Dynamical properties of a vibrating molecular quantum dot in a Josephson junction 
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We investigate dynamical transport aspects of a combined nanomechanical-superconducting de- 
vice in which Cooper pair tunneling interfere with the mechanical motion of a vibrating molecular 
quantum dot embedded in a Josephson junction. Six different regimes for the tunneling dynamics 
are identified with respect to the electron level and the charging energy in the quantum dot. In five 
of those regimes new time-scales are introduced which are associated with the energies of the single 
electron transitions within the quantum dot, while there is one regime where the internal properties 
of the quantum dot are static. 
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I. INTRODUCTION 

Inelastic scattering processes carry dynamical degrees 
of freedom which have a large influence on the electron 
dynamics. Signatures in the conductance of molecu- 
lar electronics devices^^ indicate coupling between elec- 
tronic and vibrational degrees of freedom, whereas spin 
inelastic scattering effects have been utilized in exper- 
imental studies of the magnetic properties of magnetic 
atoms and clusters f^^l A huge effort has been spent 
on studies of the influence of spin and vibrational de- 
grees of freedom on the transport through quantum dots 
(QDs)P^HII] Studies of various aspects of inelastic scat- 
tering effects are of fundamental importance. 

Recent developments towards incorporating super- 
conducting electronics into nan oelect romechanical de- 
vices open possibilities of cooling^^ISSl anc j Cooper pair 
slnittlingP2E2 In this paper we focus on the dynamical 
aspects of a QD embedded in a Josephson junction, which 
to our knowledge not has been studied previously. 

In this paper, we consider the influence of the mechan- 
ical motion of a molecular QD, embedded in a Josephson 
junction, on the supercurrent flowing across the junc- 
tion. The mechanical motion of the QD couples to the 
electronic degrees of freedom of the tunneling electrons, 
which dramatically influence the dynamics of the elec- 
tronic occupation in the QD. We can clearly distinguish 
between six different regimes, with respect to the electron 
level £o an d charging energy U in the QD, see Fig. [T] in 
which the dynamical aspects of the transport properties 
are different. In regimes II — VI, the dynamics of the 
QD properties generate new time-scales to the transport 
which are intimately associated with the energies of the 
single electron transitions in the QD. In regime I, on the 
other hand, the internal properties of the QD are static, 
even for finite bias voltages. The motion and Josephson 
current are, thus, set solely by the applied bias voltage 
and phase difference between the electrodes. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the model of the Josephson junction in which 
the molecular QD is embedded, and we derive the basic 



expressions for the Josephson current, the two-electron 
tunneling process, and the occupation numbers in the 
QD. The results are analyzed in Sec. |III| and we con- 
clude the paper in Sec. |IV| 



II. MODEL OF THE JOSEPHSON CURRENT 

We start by considering a molecular QD embedded in 
a Josephson junction, where the QD is exposed to me- 
chanical oscillations which are modeled with a Hook's 
law constant k c , and we assume that there is a bias volt- 
age applied across the junction. We illustrate the system 
schematically in Fig. [2] 

The Hamiltonian for the set-up is expressed by 



H=H l +Hr + Hq D +Ht, 



(1) 



where^ x = Ek CT ex e kcL c k<7 + EkGx( Ac L l '-ks 
X = L, R, are usual s-wave BCS Hamiltonians, whereas 

Hqd = J2a £ <ydi-d a -\-Un^n^ defines the QD, with single- 
electron levels e a = £o + ctB/2 which are spin split by 



-H.c.) 




FIG. 1: Phase diagram of the dynamics of the molecular QD, 
embedded in a Josephson junction, with respect to the energy 
level £o and charging energy U. Here, A pq = E p — E q and 
(E ,E 1 -,Ei,E2) = (0,et,e4.,e-i- + + U), where s a = e + 
oB/2, with B — 0, whereas |A| is the superconducting gap 
in the electrodes. 
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FIG. 2: Schematic view of the mechanically and electronically 
coupled QD to the superconducting electrodes (SC L) and (SC 
R). The QD is suspended on a cantilever which is modeled as 
a harmonic oscillator with spring constant k c and mass m c . 
The device is biased with the voltage V. 

the effective field B, and with charging energy U. Fi- 
nally, Ht = J2ka ^x( c Lr^ CT +^- c -) accoun ts for the single- 
electron tunneling between the lead x an( i the QD with 
rate t x . The local vibrational mode of the island is in the 
linear coupling regime given by 

t L =t { °\l + a L q), t R =t<g\l+a R q), (2) 

where &l(r) describes the coupling between the tunnel- 
ing electrons and the vibrational mode corresponding the 
to left (right) tunnel junction. The quantity q is the 
displacement operator for the oscillator. The tunneling 
matrix element t x is exponential in the displacement q, 
thus, the assumed linear coupling is a good approxima- 
tion for small q. This allows evaluation of in terms 
of the tunneling matrix elements and their distance de- 
pendence. We assume here a very general equilibrium 
geometry with no particular symmetry being required. 
The equilibrium position (q = 0) for the QD within the 
junction may be placed anywhere in between the leads. 

A. Derivation of the current 

We derive the basic formula for the supercurrent 
through the system using an analogous procedure as in 
Refs. 23 24. The total current flowing between the QD 
and lead \ can be fundamentally expressed as I x ot ^ (i) — 
- e9 t(J2 P ae x c L c k<r): from which we obtain 

4 to *)=2eRe f ({\A x {t),A^)])e^*-^ 

J— CO \ 

+ ([A x (t),A x (t')})e^ t+t '^dt'. (3) 
Here, we have defined the current operator 

M*) = E *x(*)cL(*K(*). ( 4 ) 



using that 4 CT (*) = e^^^c^e"^*-' 1 *^* and 
analogously for the ci-operators, where p, x and N x — 

Skerey c kcr Cko ' * s the chemical potential and total num- 
ber of electrons in lead x, whereas we put the chemical 
potential for the QD to zero, i.e. [Iqd = 0. 

The total current comprise one component of single- 
particle, or electron, current which is contained in the 
first term of the above expression, while the second term 
describes the supercurrent. In the present paper, we are 
interested in the supercurrent only and, therefore, we 
discard the first contribution to the total current in the 
following discussion. 

The average ([A x (t),A x (t')]) is decoupled into the av- 
erages e.g. (cj t(7 (t)cj c ^(t')) and (d a (t)d s (t')} of which the 
former is further handled via e.g. Bogoliubov-Valatin 
transformation c ko - = u k 7 k(T — ^kTW' suc h that (t = 
i-f) 

4£(mo =hx4(*)4(*')> 

=*™ k M/(-£k)e- jBkT - f(E k )e iE ^}, 

(5a) 

=i*uiv k [f(-E k )e iE * T - J(£ k )e- lBkT ]. 

(5b) 

Here, E k = \/ (e k — A^x) 2 + l^xl 2 defines the quasi- 
particle energies, whereas u k = \/(l + [e k — p, x ]/E\ s )/2 
and u k = \/(l — [e k — ^ x ]/i? k )/2, which satisfy |u k | 2 + 
|i> k | 2 = 1 and u^v k = \A X \e i ^/(2E^), with the macro- 
scopic phase <p x . We shall proceed at low tempera- 
tures, such that we can approximate /(-E k ) ~ and 
f(-E k ) « 1. 

The average (d a {t)do-(t')) is handled by transforming 
the QD Hamiltonian into diagonal form using the expan- 
sion d a = X 0a + a X s2 for the operators X™ = \p)(q\, 
A P9 = E p - E q , p,q e {0,cr,2}. 25 , such that Uqd = 
Ep=o, CT ,2 E p X pp , (E , Ef, E±, E 2 ) = (0,e^,e^,e^ + + 
U). Then, since d a (t)d a {t') = (X 0a + aX s2 )(t)(X 0rf + 
<jX a2 ){t') 1 we define the anomalous averages for the QD 
according to 

F> a (t',t) =(-i){d s (t')d a (t)) = (-i)aN 02 e iA *° r , (6a) 
F<(t', t) =i{d s {t')d a {t)) = iaN Q2 e iA ^ T , (6b) 

where 7V 02 = (X 0a X a2 ) = (X 02 ) is the average rate for 
the two-electron transition X 02 = |0)(2| at the energy 
A 2 o = E 2 — So- 
Using the prescribed procedure, we find that the su- 
percurrent between the superconducting lead x an d the 
QD at low temperatures can be written 

I x (t)=-Re J N 02 [J x { f i L )(l + a x q) 2 sin(w x i + <j> x ) 

- r x(M x )(l + a x q)a x qcos(u; x t + (f> x )]e~ lUT dt' ^ , 

(7) 
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where lo x = 2/i x , and where we have assumed the local 
approximation t x {t') ~ t x (t) — rt x (i), which is justified 
since the vibrational motion is much slower than the elec- 
tronic tunneling processes. The amplitudes J x and T x 
of the Josephson current in absence and presence of the 
coupling to the vibrational mode are given by 

J x ( Px )=2e X 2 ]g + A g0 ) 

-£ k+ (u + A 2CT )], (8a) 
IA ll/ (0) l 2 

r x ( Mx ) =2e£ ™ 1 [£L(o; + A,o) 

-4 + (u; + A 2CT )], (8b) 

respectively, where £k±(^) = l/( w — I^x ^ -^k]), k G 

The quadratic dependence of the displacement g in the 
current, Eq. Q, is justified since the q 2 <^ q for small q. 
The quadratic component, thus, merely provides a minor 
modification to the linear displacement. 



B. Two-electron tunneling process 

The dynamics of the average N 02 = (X 02 ) is calculated 
through the equation of motion 

S2 \r, -\o~ i Vx t 



{id t ~ A 2a )N 02 = -J2 (jt x((. x °' T + <?X s2 )c ks )e 



(9) 



where the correlation function is treated by using pertur- 
bation theory, which to first order in t x gives 

{[{X a ° + aX s2 )c ks ]{t)) = 

= i f t x ([G> a (t,t') + G> 2 {t,t>)]F< aS {t> ,t) 

J — CO 

- [G< ff (t,0 + G5 2 ( t ,t')]J-^(t',t))e-^ t '^. 

(10) 

The lesser and greater GFs for the QD, are here given by 



G< a (t,t') = iN a e- lA °° T , G> a (t,t>) = -iN e 



(11a) 

G< 2 (t,t>) = iN 2 e~ iA ™ r , G> 2 (t,t') = -iN s er^ T , 

(lib) 

where N p — (X pp ), denotes the occupation number for 
the QD state \p). whereas the anomalous GFs for the 
leads are expressed as 

=zau k v k [f(E k )e^ - f(-E k )e~ iE * r ], 
O =i{c ks it')c ka {t)) 

=i<7u k v* k {f(E k )e- lE * T - f(-E k y E * T }. 



The equation for Nq 2 is, thus, given by 
(id t -A 20 )N 02 = -J2 A(l + a*ff)[(l + a x q)U x (u, t') 



ia x qV x (w, t')}^- 1 ^- 1 ) —dt'e~ l{ ^ t+ ^\ 



2ir 



where the amplitudes U x and V x are 



(12) 



U x (w,t)=- 



\AJ\t 



(0)|2 



kt?£X 



2E k 



N s £ k _(A 2s -u) + N a C k+ (A a0 - u) 
■ N 2 C k+ {A 2s - lo)] (13a) 



V x (w,t)= 



\AJ\t 



xlrx 



(0)|2 



2 Ei. 



-[N £l_(A a0 -uj) 



fecrGx 

+ N s C 2 k _(A 2s -u) + N a C 2 k+ {A aQ - to) 
+ N 2 C 2 k+ {A 2s - uj)]. (13b) 

Those amplitudes describe the cotunneling processes in 
which two electrons are either added or removed from 
the QD in absence and presence of the coupling to the 
vibrational mode, respectively. 

C. Quantum dot occupation numbers 

It is clear that the supercurrent is to a great extent 
determined by the time-evolution of the electron occu- 
pation N p in the QD. We obtain those from the den- 
sity matrix p(t) = {(X pp ')(t)} pp < requiring J2 p n p = 1 - 
Within the employed level of approximation, p is de- 
termined from the master equation <9 t N = J2 X ^'"-'xN, 
where N = (N A^ N± N 2 ) T , whereas 



r x 

_A* _p x -i- A x 



iv 04. 







r x 
1 04. 








r x 

1 1-2 
pX 

X 



and where 



8(A X - 

r* =2^ACA X — ^— gg— 
pq x qP x/{A^p 



|A X I 2 

A x =2^A/ X A X — r ~ qp |Ay 'V 



tf(-A x p 



(14) 

(15a) 
(15b) 



Here, Af x is the density of electron states in lead \-> 
whereas A x ? = A pq — uj x . 



III. DYNAMICS OF THE QUANTUM DOT 

It can be seen in Tabs. |T]and|lTJ and Fig. [I] that there 
are several regimes in the (eg, C/)-space, in which the 
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TABLE I: The three regimes in which at least one transition TABLE II: The three regimes in which both transition ener- 
energy lies within the superconducting gap. gies lie outside the superconducting gap. 
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time-dependence of the quantum dot occupation num- 
bers is very different. Here, we analyze each of the 
regimes in order to elucidate the characteristic time- 
scales involved in the electro-mechanical dynamics of the 
QD. 



A. Regime I: Both transitions within the gap 

Beginning with regime I, we find that the occupa- 
tion numbers are constants of motion since the matrix 
U x = 0, hence, 9tN = 0. The non-dynamical charge 
distribution on the QD, implies that the QD motion can 
be linearly excited by varying the bias voltage and su- 
perconducting phases. Integrating Eq. (12), we obtain 



^02 =J2{ [t 1 + a xl) 2 U' x - a x {(l + a x q)q + a x q\V' x ' 
x 



cos(w x i + (j> x ) - a x q (1 + a x q)V' 



+ 2(1 + a x q)U' x ' - 2a x qV"' sin(uj x t + </> x )j 



(16) 



where U' x 



U X /(A 20 - u x ), V" = U X /(A 20 - w x f, and 

U x = U x/( A 20 - w x) 3 > and analogously for V x , V£, and 

V'" 

x 

An order of magnitude estimate for the relative ratios 
between the parameters Jg, Tg, U Xl and V x gives 



eV 



h J v 



eV 



(17) 



We note several things in this formula for the Joseph- 
son current. The first term would arise when the QD 
is rigid in space since there is no coupling to the vibra- 
tional motion. All other terms, which are proportional to 
a x , arise due to the movability of the QD in the Joseph- 
son junction. There are two terms directly depending 
on the pos ition q of the QD, as expected from previous 
studiesP222 More interesting, however, is the dependence 
on the velocity q and acceleration q of the QD, displayed 
in the fourth, fifth, and sixth terms. These dependences, 
hence, open possibilities for novel measurements and ap- 
plications of shuttling QDs. 



Following the procedure introduced in Refs. 23 24, we 
derive the effective Hamiltonian Hj by requiring that 



2edWj /d<j) x = I x , where I x is given in Eq. (18). We 



briefly discuss the salient features of the result, while a 
detailed derivation will be presented elsewhere. Within 
the given level of approximation, the Hamiltonian is lin- 
ear in the coordinates q, q, q. Hence, the classical equa- 
tion of motion from the Hamiltonian Tl osc — p 2 /(2m) + 
k c x 2 /2 + Hj is expected to be simple, in the sense that 
all features derived from the superconducting current act 
as a driving force F(t). The equation of the QD mo- 
tion will, therefore, be that of a driven oscillator, i.e. 
m c q + 7tv9 + k c q = F(t), where the damping factor 7^ 
contains the external damping due to mechanical friction. 
We notice en passant that a quadratic approximation in 
a x provides contributions that add a Josephson stiffness 
and damping. The linear spatial motion of the QD is, 
however, given by 



As we are considering small bias voltages, eV <§C |A X |, 
we proceed by neglecting the T x terms, and since we are 
considering weak coupling between the vibrational mode 
and the electronic degrees of freedom, we only keep terms 
at most linear in a x . The current, thus, reduces to 



hi*) =Jx l^G 1 + 2« x9 + 2a x ,q}U' x , - a x ,qV^,) 
x' 

x cos(w x 'i + <j> x >) + a x >q[2U x , 



x sin(u; X 'i + 4>x')W sin(o; x i + X ), 



(18) 



q{t) =q sin(£ t + S )e~'<^ 2m - q £ ^{^K 

x x ^ 

- 6 ^[3U x -2u; x (2U x , -V x \)]H(u; x , ( t> x ) 

+ J2 ( 2 xx'K')^K'^x') - ^Q XX 'K' -^x) 

x'#x 



x H{u x , -u) x ,<j> x > - 4> x ) - T.Qxx'i^x' + w x) 



xH{u x/ +u x ,cj> xl +(j> x ))}, (19) 
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where 

H{w,4>) 
whereas q v = 



■2a x ,]U^,+a x ,uj[2U^-V^] 
1 - (w/o; ) 2 



[i - (w/w ) 2 ] 2 + bW^c] 2 

7jvw sin(w£ + 0) 



cos(a;t + . 



fee 1 - (w/w ) 2 



q a x 



and A" v 



v x 

less parameters, whereas Ej 
energy. The first term in Eq. ( 19 ) describes the unper- 



k c q 2 /E* arc dimension- 
_Jy/2e is the Josephson 



turbed damped motion of the QD around its equilibrium 
position q cq , where qo and 5q are to be determined from 
the initial conditions, whereas ljq = yu% — [7jv/2m] 2 
(ujq = \/k c /m c ) is the eigenfrequency of the damped (un- 
damped) mechanical oscillations. The only time-scales 
introduced in the Josephson current are the ones associ- 
ated with the eigenfrequency loq and the Josephson fre- 
quencies lj x , which is contrasted by the properties of the 
other regimes, as we shall see below. 

This extended analysis of the expected Josephson cur- 
rent in term of the QD motion is possible only since the 
occupation numbers are constants of motion. Hence, for 
the remaining regimes we will be content with establish- 
ing the characteristic time-scales. 



Regimes II and III: Single transition within the 
gap 

Within regime II, the parameters T* Q = 0, for all tran- 

A x |, whereas 



pq 

since Af CT > 



sitions pq, and also A* 2 

Aq ct ji since A* < — |A X |. Hence, the number N 2 
is a constant of motion, while the remaining system of 
equations has the solution 



A 0<7 ft *x (*')<*«' 



(20a) 



N„(t) = f tl{t')N a {t')dt\ (20b) 



for some initial time to at which Nq = No(to). Under the 
local approximation of the tunneling rate, Nq acquires 
the time- evolution 



n 



3 A^(t<, 0) ) 3 o ;x <j{[l+a x 9][l-a x «j(i-to)]+a''j 2 (t-to) 2 /3}(i-to) 



The exponent periodically changes sign with the QD ve- 
locity q, and position q, which provides an oscillatory 
behavior of the occupation numbers No as well as for 
N a . Physically, this means that the electron density in 
the QD periodically grows and wanes as the QD moves, 
which opens for possibilities to load and unload electron 
density on the QD at different leads, or, single electron 
shuttling between the superconducting leads via the QD. 



The properties of the QD in regime III are analogously 
obtained by interchanging the roles of Nq and N 2 . 

The time-scale r for loading (unloading) electron den- 
sity on the QD is related to the transition energy A* 
( A 2(T ) , and it can be noted that r — > as A* — > — |A X | 
(A 2o . — > |A X |), while it is determined solely by the den- 
sity of electron states, r ~ 1/A/^, in the leads for transi- 
tion energies far below (above) the superconducting gap. 
The dynamics of the occupation numbers introduce an 
additional time-scale to the ac Josephson current which 
is associated with the energies of the QD transitions, and 
which is different from the one introduced through the 
Josephson frequency. 



Regime IV and VI: Both transitions below or 
above the gap 



In regime IV, all transition energies A pq < — |A 



*xl> 



which implies that T* q = 0, whereas A* q ^ 0. We can, 
then, integrate the occupation of the empty state accord- 
ing to Eq. (20a I, whereas the occupations for the one- 



and two-electron states are given by 

N„(t) = [] f eK2 1: ' tl(t " )dt " E & (f)AS>o(f)< 
x Jt " x' 

(21a) 
(21b) 



JV 2 (i)=TA ri f N a (t')dt'. 

Jtn 



While both the empty and one-electron states, again, de- 
pend on the periodic motion and velocity of the QD, the 
two-electron state also acquires an oscillatory behavior 
since it depends on the integrated time-evolution of the 
occupation of the other states. It is noticeable, however, 
that the time-evolution of the electron occupation in the 
QD here is related to the (four) time-scales associated 
with the transition energies A ct0 and A 2rT . In particular, 
the one-electron occupations N a strongly depend on the 
rate of the transitions X 0a and X a2 . The properties of 
regime VI are obtained by noticing that all T pq ^ and 
all A pq = such that the roles of No and N 2 become 
interchanged. 



D. Regime V: Gap between the lower and upper 
transitions 



In regime V, finally, A* < — |A 



and A* CT 



> IA, 



which leads to that AX < 0, T* = 0, A* 2 = 0, and 



r x 



> 0. Ao and AT 2 are, thus, determined by Eq. (20a I, 
requiring the replacements Aq^ — > r* 2 and N n —> N 2 
in the expression for N 2 . It is understood that Ao de- 
creases (increases) while N 2 increases (decreases), which 
is expected from conservation of probability. The one- 
electron occupations N a depend on the integrated time- 
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evolution of Nq and N2, that is, 

M*) = ~ E f t 2 x (t')i^N (t') - n 2 N 2 (t')]dt'. 
x Jt ° 

(22) 

The occupation of the one-electron states can, thus, be 
viewed as resulting from the imbalance between the oc- 
cupation in the empty and two-electron states 

IV. SUMMARY AND CONCLUSIONS 

We have studied the dynamics of a single level molec- 
ular QD embedded in a Josephson junction, in which the 
mechanical motion of the QD is couples to the supercur- 
rcnt. In the static limit, i.e. when the QD is rigid, the 
supercurrent between the superconducting electrodes is 
mediated via two-electron transitions in the QD, a behav- 
ior which remains also when the QD moves between the 
electrodes. The rate of the two-electron transitions nat- 
urally depends on the occupation of the QD states and 
the energies of the transitions between the states. The 
rate does, in addition, explicitly depend on the motion 
of the QD. 

Our main focus has, in the present paper, been devoted 
to extract the time-scales that are involved in the dynam- 
ics of the QD occupation numbers. It turns out that we 
can clearly distinguish between six separate regimes in 
the phase space of the single electron energy level £0 an d 
the QD charging energy U, see the phase diagram in Fig. 
[I] The boundaries between the various regimes are set by 
pairing potential (gap function) of the electrodes, and the 
dynamics of the QD occupation depends on whether the 
single-electron transitions in the QD lie within or outside 
the gap. 

In all regimes but one, there are more than two time- 
scales involved in the dynamics of the QD occupation. 



Those additional time-scales are typically defined by the 
energies of the single-electron transitions relative to the 
superconducting gap (and the chemical potential). For 
transition energies very far below (above) the gap, the 
time-scale for occupying, or deoccupying, the correspond- 
ing states are set by the density of electron states in the 
electrodes, while the time-scale tends to zero as the tran- 
sition energies approach the edge of the gap from below 
(above). As the transition energies lie within the gap 
function, the occupation numbers of the QD are con- 
stants of the motion which implies that the dynamics of 
the QD occupation depends solely on i) the two basic 
time-scales set by the bias voltage (Josephson frequency) 
and the eigenfrequency ujq of the oscillator, and ii) the 
phase difference between the electrodes. 

An experimental set-up to test our predictions would 
be possible by having gap functions | A| ~ 10 meV, which 
is pertinent to MgB2p3 in order to obtain a sufficiently 
small mechanical damping. By also acquiring a vibra- 
tional frequency of ojq/2tt ~ 1 GHz one should be able 
to tune the Josephson frequency such that ujq/loj > 0.1. 
A coupling strength oP^a/wo ~ 1CP 1 — 10 -3 should be 
sufficient for an efficient read-out. Finally, the charging 
energy U of the QD may be of the order of 1 — 10 meV, 
while the level spacing should preferably be larger than 
U in order to control the electron occupation. For this 
requirement on the charging energy, one would have ac- 
cess to the regimes I — IV and VI, whereas the regime V 
would become accessible by requiring a charging energy 
U>2\A\. 
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